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ABSTRACT 

A  method  has  been  developed  to  calculate  the  blade  frequency  point 
pressure  on  a  hull  surface  in  the  presence  of  a  propeller  and  a  free  water 
surface.  The  approach  is  based  on  that  previously  employed  in  calculating 
the  propeller-induced  hull  forces^  in  which  the  hull  is  represented  by  a 
source  distribution  smeared  over  the  hull  surface  whose  strength  is  de¬ 
termined  through  the  solution  of  a  Fredholm  integral  equation  of  the  2nd 
kind  in  terms  of  specified  propeller-induced  velocities.  Point  pressures 
are  determined  by  superposing  the  free-space  pressure  emanating  from  the 
propeller  through  its  fanning  action  on  the  pressure  due  to  the  diffraction 
potential.  A  series  of  calculations  was  performed  for  the  case  of  a  model 
tanker  tested  extensively  at  the  Norwegian  Experimental  Tank  (Ship  and  Ocean 
Laboratory,  Norwegian  Hydrodynamic  Laboratories,  Marine  Technology  Center). 
The  agreement  between  experimental  and  predicted  values  varies  from  point  to 
point;  in  general,  the  maximum  deviation  of  amplitude  is  between  ±30%,  but 
the  phases  are  in  poor  agreement.  A  new  method  is  suggested  in  which  the 
hull  is  represented  by  a  doublet  distribution  rather  than  a  source  distribu¬ 
tion.  This  approach  which  reduces  computing  time  considerably  permits 
improved  accuracy  and  results  in  much  better  correlation  between  experimental 
and  calculated  values. 
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NOMENCLATURE 


c(p)  expanded  propeller  chord  at  radius  p 
d  distance  of  propeller  axis  from  free  surface 

lm(x)  modified  Bessel  function  of  first  kind  of  order  m 

K^x)  modified  Bessel  function  of  second  kind 

^m,n^  modified  space  function  (after  chordwise  integration) 
variable  of  integration 

^,n^(p)  =  L,  -(p)  =  spanwise  blade  loading,  lb/ft 
a  9  n 

,  _  loading  coefficients 

\,n  a 

i nteger  multiple 
integer 

number  of  blades 


order  of  chordwise  modes 
pressure,  lb/ft2 

order  of  harmonic  of  velocity  field 


Descartes  distance 


propeller  radius,  ft 
surface 


x,  r.cp 


time,  sec 

free-stream  velocity 

axial  component  of  velocity  (hull  axes  system) 
variable  of  integration 

axial,  radial  and  tangential  components  of  velocity  field 
of  propeller 


nim 


.1 
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Vx,y,z 

WZ 

X ,  Y,Z 


x,y,z 


X,  r,cp 
®  =  -Ot 

eb 

A^(x) 

X 

X 


§,P,eo 

5, VC 

Pf 

a 

CT1 

To(p) 

f  ,cp 


Q 


velocity  components  in  Cartesian  coordinate  system 

vertical  component  of  velocity  (hull  axes  system) 

hull  Cartesian  coordinate  system 

propeller  Cartesian  coordinates  of  field  point 

propeller  cylindrical  coordinates  of  field  point 

=  blade  angular  position 

projected  propeller  semichord,  radians 

defined  in  Appendix 

order  of  harmonic  of  loading 

variable  of  integration 

cylindrical  coordinates  of  point  on  propeller  blade 

hull  Cartesian  coordinates 

mass  density  of  fluid,  slugs/cu  ft 

source  strength 

angular  measure  of  skewness 

maximum  thickness  of  blade 

velocity  potential 

angular  velocity  of  propeller 
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INTRODUCTION 

Propeller-induced  ship  hull  vibration  continues  to  be  a  major  uncer¬ 
tainty  to  the  naval  architect  designing  any  class  of  ship.  In  all  cases,  the 
amelioration  of  hull  vibration  is  desired,  with  different  levels  of  tolerance 
since  vibration  causes  discomfort  to  passengers  and  crew,  produces  material 
strain  and  fatigue  and  generates  noise  which  impedes  communications.  With 
today's  trend  towards  higher  speed  and  larger  ship  dimensions,  requiring 
greater  power  to  be  delivered  to  the  propeller,  the  problem  of  propeller- 
induced  hull  vibration  has  become  of  paramount  importance,  necessitating  im¬ 
provement  of  the  methods  for  predicting  propeller-excited  forces  and  moments. 

The  methods  have  to  be  reliable  and  practical  for  use  at  the  design  stage.  It 
is  believed  that  only  upon  the  development  and  use  of  a  readily  applicable  set 
of  computer  programs  can  future  vessels  be  designed  which  have  reduced  propeller¬ 
generated  vibration. 

(1  2)* 

Investigations  ’  '  have  been  initiated  to  provide  a  satisfactory  proce¬ 
dure  and  furnish  adequate  information  to  allow  the  naval  architect  to  determine 
propeller-induced  hull  excitation.  A  generalization  of  the  earlier  procedures 

for  prediction  of  propeller-generated  vibratory  hull  forces  on  arbitrary  hulls 

(3) 

has  recently  been  published.  '  The  task  remaining  is  to  demonstrate  the  reli- 
ability  of  this  method  by  comparison  of  calculated  pressures  with  E.  Huse'sv 
extensive  measurements  of  point  pressures  on  a  tanker  model  with  CB  =  0.80. 

D 

The  method  of  calculating  the  blade-frequency  point  pressure  on  the  hull 
and  the  corresponding  hydrodynamic  forces  and  moments  requires  assuming  a 
source  distribution  over  the  hull  surface  to  simulate  the  ship  in  both  uniform 
axial  and  propeller-generated  flows.  The  strengths  of  the  singularities  are 
determined  through  the  solution  of  a  Fredholm  integral  equation  of  the  second 
kind  knowing  the  propeller-induced  velocity  on  the  hull.  To  carry  out  the 
numerical  solution  of  the  integral  equation,  the  hull  surface  is  subdivided  into 
quadrilaterals.  The  quadrilateral  size  is  crucial  in  calculating  the  hydrody¬ 
namic  forces  and  point  pressures,  and  the  large  number  of  quadrilaterals 
required  to  maintain  the  desired  accuracy  together  with  the  complication  arising 
from  the  calculations  of  the  propel ler- induced  velocity,  makes  this  numerical 

Superior  numbers  in  text  matter  refer  to  similarly  numbered  references  listed 
at  the  end  of  this  report. 
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approach  expensive  and  hence  not  attractive  for  practical  application. 

As  mentioned  in  Reference  3,  a  new  approach  has  been  developed  with  MARAD 
support,  in  which  the  hull  is  represented  by  a  distribution  of  doublets  whose 
strength  will  be  determined  by  the  same  integral  equation  as  before,  but  with 
the  known  function  replaced  by  the  negative  of  the  propeller  velocity  potential 
function.  The  latter  has  no  directivity  and,  hence,  there  is  only  a  single 
value  at  each  point  on  the  hull,  whereas  the  present  method  requires  calcula¬ 
tion  of  three  components  of  the  propeller-induced  velocity  field  to  determine 
the  normal  velocity  at  each  point  on  the  hull  surface.  The  computer  time  for 
evaluating  the  strength  of  the  doublet  distribution  is  considerably  less  than 
that  currently  required  for  the  determination  of  the  hul 1 -source-dens i ty  distri¬ 
bution  and  as  has  been  seen,  greater  accuracy  can  be  achieved  within  a  shorter 
time.  This  will  be  further  demonstrated  in  a  future  study  correlating  predicted 

and  experimental  values  for  forces  on  the  Series  60  hull  and  point  pressures  on 
(4) 

Huse's  tanker  model. 

The  present  study  is  sponsored  by  the  Naval  Sea  Systems  Command,  General 
Hydromechanics  Research  Program,  under  Contract  N00014-79-C-0240,  administered 
by  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center. 
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PROFELLER-INDUCED  VIBRATORY  PRESSURES  AND  FORCES  ON  A  HULL 
A.  VELOCITY  POTENTIAL 

Let  the  velocity  potential  of  a  ship  stationary  in  uniform  flow  of  veloc¬ 
ity  U  along  the  +X  direction  (see  Figure  1)  in  the  presence  of  an  operating 
propeller,  with  and  without  free  surface,  be  denoted  by  9  .  When  the  pro¬ 
peller-induced  pressure  and  velocity  fields  are  calculated,  the  fact  that  the 
propeller  operates  in  a  spatially  varying  inflow  field  must  be  taken  into 
account.  The  presence  or  absence  of  the  free  surface  should  be  taken  into 
account  by  a  proper  image  system  reflecting  the  fact  that  the  free  surface 
behaves  either  like  a  "rigid  wall"  or  a  surface  of  constant  pressure. 

The  ship  velocity  potential  can  be  expressed  as 

9=9+9  (1) 

where  9gs  is  the  velocity  potential  of  a  ship  in  a  uniform  stream  (+U)  and 
m  o  the  potential  of  the  ship  in  the  presence  of  the  propeller-induced 
velocity  field. 

When  the  ship  hull  is  represented  by  a  source  distribution  over  the  hull 
surface,  the  first  potential  will  be  given  as 


s  s 


(2) 


where  oq  is  the  source  density  which  creates  the  hull  in  the  presence  of  a 
uniform  stream  U;  R  and  Rj  are  the  Descartes  distances  to  any  field  point 
P(X,Y,Z)  from  any  hull  point  Q(?,T1,C)  and  its  corresponding  positive  image 
in  the  free  surface;  S  is  the  wetted  hull  surface.  Taking  into  account  the 
symmetry  of  the  flow  with  respect  to  the  waterplane,  and  the  symmetry  of  the 
hull  with  respect  to  the  centerplane,  the  image  system  is  constructed  so 
that  the  "zero"  Froude  number  or  "rigid  wall"  condition  is  satisfied  under 
the  assumption  that  the  ship  speed  is  low. 


When  the  ship  hull  is  represented  by  a  source  distribution  of  density  a 
over  the  wetted  hull  surface  to  create  the  hull  in  the  presence  of  the  pro¬ 
peller  and  its  negative  image  in  the  waterplane,  then  the  ship  velocity 
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potential  is  given  as 

%  =  ‘  h  S  a(i  -  (3) 

where  R  and  R.  are  the  Descartes  distances  to  a  field  point  P  from  a  hull 

i 

point  Q  and  its  image  point  Q.;cr  is  the  blade-frequency  dependent  source 
dens i ty . 

It  is  clear  that  for  all  points  on  the  free  surface  where  R=R.  the 
high  Froude  number  limiting  condition  cpgp  =  0  is  satisfied. 

B.  INTEGRAL  EQUATION  FOR  THE  HULL  SOURCE  DISTRIBUTION 

It  is  assumed  that  the  ship  hull  is  represented  by  a  distribution  of 
sources  over  the  hull  surface  of  unknown  density  a,  which  is  determined  by 
the  requirement  that  their  normal  velocity  balance  any  incident  flow.  This 
required  kinematic  boundary  condition  yields  the  following  integral  equation 
for  the  unknown  source  density  distribution,  a(P) , 


where  represents  either  the  kernel  of  Eq.(2)  or  Eq.(3).  In  either  case, 
d  ° 

is  the  derivative  normal  to  the  hull  surface  at  point  P;  S  is  the  wetted 
surface  of  the  hull.  This  equation  is  a  Fredholm  integral  equation  of  the 
second  kind  for  the  unknown  source  densities  o,  which  are  determined  in 
terms  of  the  known  onset  flow  on  the  hull  surface  specified  on  the  right-hand 
side  of  the  equation.  The  following  flow  field  imposed  on  the  locus  of  the 
hull  is  considered:  a  flow  imposed  on  the  hull  due  to  uniform  longitudinal 
flow  and  two  components,  the  real  and  imaginary  parts  of  the  blade  frequency 
velocity  field  induced  by  the  action  of  the  propeller.  For  each  given  onset 
flow,  the  source  densities  are  computed  by  solving  Eq.(4). 

The  coded  numerical  procedure  developed  by  Hess-Smith^  is  applied  to 
solve  this  integral  equation.  The  three  onset  flows  considered  produce 
three  densities,  namely,  a  ,  due  to  the  uniform  steady  flow  and  and  o^ 
due  to  the  real  and  imaginary  parts  of  the  onset  flow  due  to  the  propeller 
and  its  image  in  the  free  surface  corresponding  to  high  frequency  Froude 
number.  The  use  of  the  negative  image  of  the  propeller  located  at  the  re¬ 
flected  position  above  the  free  surface  satisfies  the  high  Froude  number 
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conditions  on  the  free  surface.  For  points  on  the  free  surface,  this  con¬ 
dition  yields  zero  perturbation  pressure  or  equivalently  doubles  the 
vertical  component  of  the  velocity  induced  by  the  direct  propeller. 

Thus  the  inversion  of  the  integral  equation  is  achieved  by  a  numerical 
procedure  implemented  by  the  Hess-Smith  program^  provided  the  right-hand 
side  is  known.  This  part  constitutes  the  major  computing  effort  of  this 
endeavor. 


C.  PROPELLER-INDUCED  VELOCITY  FIELD 

1)  Due  to  Blade  Loading 

Considerable  effort  has  been  devoted  to  the  development  of  an  unsteady 
lifting  surface  theory  for  propellers  operating  in  a  spatially  nonuniform 
hull  wake.  The  theoretical  analysis  and  numerical  techniques  developed  at 
Davidson  Laboratory  are  based  on  modeling  the  propeller  blades  with  sources 
and  pressure  dipole  singularities.  The  computer  program  permits  the  calcula¬ 
tion  of  steady  and  unsteady  blade  pressure  distributions  and  resulting  bearing 
forces  and  moments  for  a  given  propeller  geometry  (i.e.,  number  of  blades, 
pitch,  chord  length,  camber,  thickness,  skew)  and  a  given  inflow  field.  De¬ 
tails  of  the  theory,  numerical  analysis,  and  computer  programs  have  been 

reported  elsewhere.  The  theory  has  been  extended  to  compute  the  field 

.  „  .  ...  (8,9)  ,  (10,11) 

point  velocities  and  pressures. 

The  steady  and  unsteady  blade  loadings  are  used  to  calculate  their  field 
point  inductions  as  responses  to  the  propeller  blades  operating  in  the  hull 
wake.  The  wake  survey  has  been  harmonically  analyzed  in  terms  of  the  shaft 
rotation  frequency,  and,  upon  summing  over  all  the  blades,  only  certain 
harmonics  contribute  to  the  net  exciting  forces  and  moments  acting  on  the 
shaft,  namely,  XN  for  the  thrust  and  torque  and  Xn±1  for  the  side-bearing 
forces  and  moments  where  X  is  any  integer  and  N  is  the  number  of  blades. 
However,  all  the  harmonics  of  loading  contribute  to  the  forces  on  an  individ¬ 
ual  blade  and  to  the  blade  rate  of  radiated  velocity  and  radiated  pressure 
fields  of  the  propeller  (e.g.,  at  points  on  the  hull  surface). 

The  velocity  potential  at  a  point  (x,r,cp)  (in  the  propeller  coordinate 
* 

system)  due  to  the  loading  of  an  N-bladed  propeller  operating  in  a  spatially 
*See  Figures  2,3  and  4. 
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varying  inflow  field  is  given^’^  by, 

4(x,r,cp;q)  =  Re  J  S  2  L^’  ^  (p)  e  ' qfit  K^m>^  dp 
p  n=l  X=o 

where  the  function  which  is  known  as  the  propagation  or  influence 

function,  is  given  for  q=AN,  n^X-i'N,  by 

!Kj  in'(a'-cP)  r  i^KiCcr'- 


lm(aUNIp)Km(aUNlr)A{n)(X6b)-  i  J 


i  k(x-o'/a) 


■  'BOK'P)K„(lklr)-A(n)[(n,-  *)«„>}  (5) 


for  p  <  r  (for  p  >  r,  p  and  r  are  interchanged  in  the  modified  Bessel  func¬ 
tions).  The  symbol  (x)  denotes  the  results  of  the  chordwise  integra¬ 

tion  which  differs  with  the  assumed  chordwise  modes  (see  Appendix). 

The  velocities  in  the  x,  y,  and  z  directions  of  the  Cartesian  coordinate 
system  with  origin  at  the  propeller  center  at  a  point  x,r,cp  due  to  the  pro¬ 
peller  operating  in  spatially  varying  inflow  will  be  given  by  the  appropriate 
components  of  the  gradient  of  the  velocity  potential  4. 

The  axial  component  in  nondimens ional  form  is 


(6) 


and  the  y  and  z  components  are  derived  from 


V  4  4  .4 

-jL  m  -X.  s  .  — L  s  i  ncp  -  —  coses 

U  U  U  r  U  Y 


4  4  ,4 

—  b  — JL  coscp  -  —  sincp 
U  U  ^  r  U 


(7) 


where  4  and  4  are  given  by  the  derivatives  with  respect  to 

r  cp  ' 

of  the  velocity  potential  (Eq.5). 


r  and  cp 
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In  nondimensional  form  in  the  steady  state,  nonuniform  inflow  case, 
q=j£*0,  m=X 


^ _  V°> 

u 


—  -  J  S  S  L<^n)(p)K<X’n>dp 

p  n=i  X=0 


4-  "-4^"  I  s  S  l-<X'")(P)K'X>f;,dp 

U  U  p  n-1  X=0  r 


(8) 


•4—  -  '  7  J  -E 

o  p  n=l  X=0  Y 


U 


and  in  the  unsteady  blade  frequency  case,  q=N,  1^1=1  and  m=X-N 


v!N)  *V(N) 


vfN)  i  (N) 


2  2  • 

M*.") 

(p)K^“N,n) 

n= 1  X=0 

+  conj  . 

(p)k<x+n>S) 

2  2  • 

rL(x»n) 

(pjjjU-N.n) 

n=l  X=0 

r 

+  conj .  | 

r< 

-» 

l  1 

,p)R(x*N,5) 

(9) 


o  p  n=l  X=0  Y 


U 


+  conj.  [L^,n^(p)K^X+N,n^]  }dp 


where  it  is  seen  that  the  propeller  induced  velocity  field  is  made  up  of  a 
large  number  of  combinations  of  the  frequency  constituents  of  the  loading 
function  with  those  of  the  so-called  "propagation  function"  |n 

uniform  inflow  the  only  blade  loadings  present  are  at  X“0  .  The  correspond¬ 
ing  velocities  are  obtained  by  restricting  Eqs.(8)  and  (9)  to  the  X«0  terms. 
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The  program  requires  that  the  field  points  be  in  cylindrical  coordi¬ 
nates  x,r,cp  with  origin  at  the  propeller  center.  The  axial  x  positive 
aft,  and  radial  r  distances  are  in  terms  of  propeller  radius.  The  angu¬ 
lar  ordinate  cp  is  measured  counterclockwise  from  the  vertical  (12M) 
position  looking  forward  for  the  direct  system.  The  connections  between 
cylindrical  and  rectangular  propeller  coordinates,  with  y  positive  to 
starboard  and  z  positive  upward  are  for  the  direct  system;  r  =  Jyz+Zs  » 
cp  =  tan  '(-y/z),  -tt  £  cp  £  tt.  For  the  image  propeller,  r.  =  »/y2+(2d-z)2  . 
cp.=  tan" '  (-y/(2d-z))  measured  clockwise  from  the  6  o'clock  position,  where 
d  =  distance  of  the  propeller  axis  from  the  free  surface.  For  the  image 
propeller  the  formula  for  Vz/U  in  Eq - ( 7)  is  opposite  in  sign. 

The  propagation  functions  were  given  first  in  References  8  and  9. 

They  were  given  later  in  Reference  11,  the  documentation  of  the  computer 
program  for  the  propeller-generated  velocity  field,  with  some  changes  in 
sign  to  agree  with  the  sign  convention  of  Reference  10  for  the  pressure 
field.  They  are  repeated  here  with  the  correct  signs^^  to  obviate  con¬ 
fusion. 

;)  For  the  axial  velocity  component  Vx 

i  ( \a 1  -mcp-aXNx) 

K  (m,n)  =  +  — - 

x  4npfu=aro 

•jaXN(a2£N  -  m/p2)  I m(a |  XN|p)Km(a|  XN|  r)A^  {XBj 

-  i-  f  ilM — ~  duj'  for  p  ^  r  (10) 

o 

where  h(u)  *  (u-aXN) (au-a2XN+  m/p2) |  ( | u-aXN| p) 

•  Km(  I  u-aiN  I  r)e '  u  u/a)6b]  . 


The  integrand  at  u=0  is  evaluated  in  Reference  8  as 
(a)  When  q*X=0,  m»X*0 


1  im 
u-0 


0 
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« 


(b)  When  q=X=0,  m=XftO 


m.(n) 


lim  Fh(u)-h(-u)  1  f(1/p  )(p/r)  A;.!(X6b) 

L  u  J  ""t(  1/p2)  (r/p)mA  n  (^6b) 


for  p  £  r 
for  p  ^  r 


(c)  When  q=iN^0 


1  im  j~h(u)  -h(-u) 

u-*  0  L  u 


-  -21  lm(a|  ^N|p)Km(a|  XNI  r)  j 


-  J  [ 2a2  j£N  -  m/p2  -  i  (x-a'/a)a£N(a2£il-m/p2)  (^©b) 

-  i6b^(a2£N-m/p2)A^  (Xeb)  i 

+  |pKm(a|  £N|r)[lm_1(a|XNlp)  +  I  m+1  (a|  £NI p)  ]  -  MjaUNIp) 

*  C  Kjp.  |  (a  I  XNI r)  +  Km+,  (a|  iNI  r)  ]  (  •  aXN(a2£N-m/p2)A^  (\6b)  ,  for  p^r 

where  the  upper  sign  (-)  is  used  when  >£>0  and  the  lower  sign  (+)  when 

X<  0.  When  p>r,  p  and  r  are  interchanged  in  the  factors  in  braces  (involv- 

(n) 

ing  modified  Bessel  functions  and  their  deri vati ves)  . For  Aj  (x) ,  see  Appendix. 
The  closed  term  of  (10)  can  be  shown  to  be  equal  to  zero  when  q=£=0  . 

ii)  For  the  radial  velocity  component  Vr 
When  p  <  r 

i  ...  i  (Xa'  -mcp-aXNx) 

r  4np,u  ar 


fu  o 

{-  ~  a  I  XN I  ( a2  iN  -  m/p2)  lm(al  XNI  p)  [«„,_!  (aUNI  r) 

,  »  g,  (u)  -  g1  (-u) 

J  u 

o 


(~\  1  ”  Si  V“/  -  Si  V--/  > 

*  Km+,(»liNlr)A(n,(Wb)  J  — - “  -  do} 


(Ha) 
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where  g  ^ ( u)  =  I u-aiNI (au-a2iN+  m/p2) lm( lu-aXNI p) 

•  [Km_1  (lu-aXNl  r)  +  Km+1  ( I  u-a£NI  r)  ] 

.  <.iu(K-a'/=)A(S)t(x.  u/a)eb] 


It  can  easily  be  shown  that  when  q=£=0  the  closed  term  of  (71 a)  becomes 

(i/2)(p/r)m(m/rp2)A(n)(Xeb) 

where  m=\  since  q=0  . 

/0\ 

The  integrand  at  u=0  is  evaluated^  'as 

(a)  When  q=£=0,  m=X=0 


1  im  T 
u-oL 


9,(u)  -  9j(-u) 


]  =  (^a/r)A(n)(0) 


(b)  When  q=£=0,  m=X/0 


11m  -  9l(u)  ~  gl('u)  1 
u-oL  u  J 


(2/r)(p/r)m 


•  j[a+(  im/p~)  (x-a»/a)  Ja^^  (X0b)  +  ( im/p2)  (0b/a)Ajn^  (X6b)  } 


(c)  When  q=£Nj*0 

;im0  [  9l(U)  u  9l('U-  ]  =  ±2hm(al^!p)[Kni.,(aUNIr) 

+Km+i(a|ANIr)JJ  j[ 2a2£N-m/p2  -  i (x-a'/a) (a£N) (a2XN-m/p2]A(n) (X9b) 
-ieb£N(a2£N-m/p2)A]^  (\0b)  j  +  a£N(a2£N-m/p2)A^  ( \6b) 

*  1  Pt  lm-i  (alXNI  p)  +  lm+1(a|£Nlp)][Km_,(a|iNlr)  +  K|n+,  (a|  £NI  r)  ] 

-  r  I  m  ( a  I  XN I  p)  L  K|n_2  (a  I  £NI  r)  +  2Km(a|£Nlr)  +  »W2(a|  XNI  r)  J } 
where  the  upper  sign  (+)  is  used  when  £>0  and  the  lower  sign  (-)  when  £<  0  . 
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When  p  >  r 


-  i  (Xc'-mcp-aXNx)  r 

K(m,nMp>r)  =  +  — - - -  •  i+(i/2)al£NI  (a2£N-m/p2)K  (al^ilp) 

r  4ttP,u  ar  1  m 

f  o 


•Cl  ,(aUNIr)  +  l  (a|iN(  r)  ]A(n)Ueb) 


+(1/2tt)  J 


g2(u)  -  92(_u) 


(lib) 


where  g2(u)  =  -|u-a£Nl (au-asXN  +  m/ps)Km(lu-a£NI p) 


.[  Im_j  (I  u-ajLNI  r)+lm+1  (lu-aJJNI  r)  ]-e' U^X  °  ^A^L  (X-u/a)6bD 


In  this  case  when  q=i=0  ,  the  closed  term  is 

-(i/2)  (r/p)m(m/rp2)A^  (X8b) 
r  o\ 

The  integrand  at  u=0  is  evaluated'  ;  as 
(a)  When  q=X=0,  m=X=0 


-  g7(u)  -  g2(-u) 

1  im  |  _£ _ £ 


(b)  When  q=j£=0,  m=Xi<0 

,  r  g,(u)  “  g7(-u)  n  . 

J: - ~_i -  j  =  (-2/ r)  ( r/p) 


*  j[a+(  im/p2)  (x-c ’/a)  Ja^  (X8.)  +  ( 1  m/  p2 )  ( ©b/ a)  A  j  ( ^®b)  j 
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(c)  When  q=£N^O 

r  g2(u)  -  g2(-u) 

i?0L - u -  J  =  +2Km(a|  2NI  p)  [  (a|  £NI  r)+lm+1  (a|  £NI  r)  ] 

•  i[2a2JJN-m/p2-i(x-a'/a)a£N(a22N-m/p2)]A(n^(Xeb) 
-i6bAN(a2£N-m/p2)A[n^(xeb) j  +  (aXN) (a2 AN-m/p2)A^n) (X0b) 

•  lp[Km_1(a|£Nlp)  +  Km+](aUNIp)][lm_,(a|£Nlr)  +  lm+J  (a|  iNI  r)  ] 

-  rKm(a|lNlp)[lm_2(a|iNlr)  +2lm(a|XNIr)  +  I  m+2  ( a  J  £NI  r)  J } 
where  the  upper  sign  (-)  is  used  when  ^>0  and  the  lower  sign  (+)  when  £<0  . 
When  p=r  ,  then 

j^(m» n)  ( p_r)  _  [  £q  .(11a)  +  Eq.  ( 1  lb)]  evaluated  at  p=  r  .  (lie) 

iii)  For  the  tangential  velocity  component  V 

-i  i  (Xcr'-mcp-aXNx)  r  /-■> 

kJT’  *  +  - ~2 - {m(a2jeN-m/p2)lm(al£Nlp)Km(aUNIr)Al  '  (X0b) 

+  {'/^J  ml_  f(u)  ~  f(~u)]du}  (12) 

o 

where 

f  (u)=(au-a2AN+m/p2)  I  m  ( I  u  -  a  XN  Ip)  Km  ( I  u  -  a  XN I  r)  -e'1^*  °  ^  A^C  ( X-u/a)  0b] 

It  can  be  shown  that  when  q=i^0  ,  the  closed  term  of  (12)  becomes 

(-m/2p8)  (p/r)mA^^  (X0b)  for  p  s  r 

(-m/2p2) (r/p)mA^n^ (X©b)  for  p  s  r 

The  integrand  of  (12)  evaluated  for  u=0  yields 
(a)  When  q«£*0,  m*X*0 

iimm|~ftu)  -  tol  1  ■  o 

u-0  L  u  J 
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(b)  When  q=£=0,  m=X>40 

lim  ml  ~  1  ■  }[a  +  i  (x-a'/a)m/p2]A^  (X6.  ) 

u-  0  L  u  J  ° 


+(ieh/a)(m/P2)Ajn)(xeb)i  f  , 


(p/r)m  for  p  s  r 


’p)m  for  p  ^  r 


} 


(c)  When  q=£Nfi0,  m=X-lN 


Tf(u)  -  f(-u) 

1  im  m  — ‘ - » — L 

L  u 


u  “*  0 


+2n'.!.:  a-,  •  o'/a)  (a22N-m/p2)A(n)  (X0b) 


-(iGb/a)  (a2iN-m/p2)Ajn^  (X0b) }  I  (a|  JtNI  p)  K  (a|  XNI  r) 
±m(a24N  -  m/p2)A^^  (X9b) 

*  i  pKm(al  XNI  r)[  lm1  (al  XNI  p)  +  1  (n+,  (a |  XN I  p)  J 

-r  lm(a|  XNI  p)  [Km_,  (a|  -fcNI  r)  +  Km+J  (a|  XNI  r)  ]  ) 


for  p  £  r  ,  otherwise  p  and  r  are  interchanged  in  the  product  of  the 
modified  Bessel  functions  and  its  derivative  (the  factor  in  braces  in  the 
second  term).  Here  the  upper  sign  (+)  is  used  when  l>  0  and  the  lower 
sign  (-)  when  i  <  0  . 

When  m=0,  i.e.,  X=£N  whatever  the  value  of  X, 


2)  Due  to  Blade  Thickness 

In  considering  the  blade  thickness,  it  must  be  remembered  that  the 
marine  propeller  with  its  blade  describing  a  helicoidal  surface  is  one  of 
the  few  nonplanar  lifting  surfaces.  The  blade  thickness  plays  a  dual  role, 
influencing  both  lifting  and  nonlifting  characteristics  of  the  blade.  As 
part  of  the  nonplanar  surface,  it  induces  a  continuous  component  of  velocity 
on  points  of  the  surface  itself,  thus  affecting  the  blade  loading  distribu¬ 
tion.  The  symmetrical  flow  disturbance  caused  by  the  blade  thickness 
influences  the  pressure  distributions  on  the  suction  and  pressure  sides  of 
the  blade  without  contributing  to  the  net  loading.  The  so-called  "nonplanar 
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thickness"  has  been  shown  to  have  small  effect  on  the  blade  loading  and 
hydrodynamic  forces  and  moments  and  presumably  negligible  effect  on  the 
velocity  and  pressure  fields  in  the  neighborhood  of  the  operating  propeller. 
The  so-called  "symmetrical  flow  disturbance"  thickness  has  an  appreciable 
effect  on  the  blade  pressure  distribution  and  a  sizable  effect  on  the 
velocity  field  as  well. 

In  this  present  study  the  effect  of  the  "nonplanar"  thickness  has  been 
neglected  in  evaluating  the  propeller  loading,  but  the  effect  of  the  sym¬ 
metrical  blade  thickness  distribution  is  included  by  using  the  thin-body 
approximation.  The  blade  cross-section  is  approximated  as  an  intersecting 

of  two  circular  arcs  (i.e.,  lenticular  cross-section)  and  it  is  represented 

(9) 

by  a  source-sink  distribution  smeared  over  the  blade  area.  '  This  approxi¬ 
mation  was  shown  to  be  justifiable  for  calculating  velocities  at  points  in 
the  field  but  not  on  the  propeller. 


The  velocity  potential  due  to  the  blade  thickness  is  thus  given  as 


j6=-® 


[ s i n ( xeb/a)  -  ( X©b/a) cos ( X0b/a) ] 


{ I  XN(  IX  +  aiNlp)K£M(lX  +  a£Nlr)e 


i  X(x-a'/a) 


-I  ^(l-X+ajWlpjK^d-X+ajW1  r)e"'X^X“CT  /a^  jdXdp  ,  for  p<  r 

03) 

The  velocity  components  along  the  x,  y,  and  z  coordinate  axes  will  be  given 
by  Eqs . (6)  and  (7). 


It* 


i)  Axial  velocity  component,  V ^ 

Taking  the  derivative  of  Eq .(13)  with  respect  to  x 

VTX  _  3$(x,  r,cp;t)  _  2aN  v  -i£NQt  iiN(ax-cp) 

U  USx  “  6 

•J  ■^-+e^p2^  J  F(X)[(X  +  aJtN)G(X.Jt) 

+  (X  -  aXN)G(-X,je)]dXdp  (14) 


where 

F(X)  =  [sin(X6b/a)  -  (X0b/a) cos (X6b/a) ]/X2 
G(X,£)=  l£N(IX  +  a^NI  p)  K^(  I X  +  a£NI  r)e' X(x_CT'/a) 


G ( — X,  £)-  lXN(l-X  +  aXNI  p) KXN( !  -X  +  a£NI  r)e“'  X(x"a'/a)  (15) 


|t  can  be  shown  that  the  integrand  is  zero  when  X=0  .  Equation  (14)  is 
written  for  p<  r  .  For  p> r  ,  p  and  r  are  interchanged  in  the  modified 
Bessel  functions. 

In  the  unsteady  blade-frequency  case,  1^1  =  1,  the  axial  velocity  is 
(cf  Reference  8) 


v  <-"> 

-  +  conj  . 


2  V 


(-N) 


TX 


U 


-1) 


4aN  -iN(ax-cp)  r>  ( 1  +  a2  p3)  tp(p) 

IT8  6  p  eb 

00 

•J  F( X) [ ( X-aN) G(X, -1)  +  (X+aN)G(-X,-l) JdXdp  (16) 


In  the  steady  state  case,  l  *  0 


=  +  iiiH  f  0  +  38 p2) 
u  tt®  jp  eb  c (p) 

00 

*  J  XF ( X)  l0(Xp)K0(Xr)cos[X(x-o'/a)]dXdp 


07) 


R-2 


i)  Tangential  velocity  component,  V, 


The  nondimensional  tangential  velocity  is  given  by 


VTcp  _  _  2aN  £  iNc-iiN[flt-(ax-cp)]  r»  ( 1  +  a2p2] 
U  "  "  ^  i— «  {  efa 


*J  F(X)[G(\,£)-G(-X,i)]d\dp 


The  corresponding  blade  frequency  component  is  given  by 


■J  F(X)[g(X,-1) -G(-\, -1) dXdp 


In  the  steady  state  case 


ii)  Radial  velocity  component,  VTr 

The  nondimensional  blade  frequency  radial  component  will  be  given  as 


•or  p  <  r; 


VTr  S§r(x,r,<p;t)  _  2aN  •  -iiNntei£N(ax-cp) 

- - - -  —  6 

P  (1  +  a2p2)  ^(p)  °p  r/^ 

.  {'JM(n^ip)aKjN('feT'N'r)  «IMx-°'/a) 
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U 


(p  =  r)  =  [Equation  (22a)  +  Equation  (22b)]  evaluated  at 


P  =  r  (22c) 


In  the  unsteady,  blade-frequency  case  |i|=l,  the  radial  component  is  obtained 
as 


Vt 


(-N) 


VT 


(N) 


(-N) 


+  conj 


=  2 


(l  =  -1) 


(23) 


that  is,  twice  Eqs.  (21a) ,  (21b) ,  or  '"’lc)  with  £=  -1. 

A  computer  program  has  been  developed  and  adapted  to  the  CDC  high-speed 

(8,9,11) 

digital  computer  for  the  evaluation  of  the  field  point  velocities, 
provided  the  propeller  blade  loadings  have  already  been  evaluated  through 
another  numerical  code  designated  as  PPEXACT^  Thus,  the  necessary  informa¬ 
tion  for  the  solution  of  the  integral  Eq.(4)  has  been  established. 


D.  SOLUTION  OF  THE  INTEGRAL  EQUATION 

Equation  (4)  is  inverted  by  using  the  Hess-Smith  program  developed  at 
HcDonnel 1 -Dougl as  aircraft  Corporation.  It  should  be  noted, however, that  the 
hull  and  its  boundary  conditions  are  referred  to  a  new  cartesian  co¬ 
ordinate  system  attached  to  the  hull,  as  shown  in  Figure  1.  The  X-axis 
is  in  the  same  direction  as  x  attached  to  the  propeller,  the  Z-axis  is 
towards  starboard  as  is  the  y-axis  in  the  propeller  coordinate  system,  and 
the  +Y-axis  is  vertically  downwards  opposite  in  direction  to  the  z-axis  of 
the  propeller  system  with  the  XZ-plane  taken  to  be  the  waterplane.  The 
relation  between  the  two  coordinate  systems  is  taken  into  consideration  in 
the  program.  Only  the  hull  coordinates  of  points  on  the  hull  surface  are 
specified  as  input  to  the  program.  Whenever  there  is  a  need  to  make  use  of 
the  corresponding  propeller  coordinates,  those  are  determined  automatically 
by  the  program. 

The  hull  surface  is  subdivided  into  a  large  number  of  quadrilateral 
elements,  each  of  which  is  considered  to  be  a  source  panel  of  uniform 
density.  The  boundary  conditions  (i.e.,  the  propeller-induced  point  veloc¬ 
ity  field)  are  specified  at  the  null  point  in  each  of  the  source  panels, 
thus  generating  as  many  linear  algebraic  equations  as  there  are  unknown  sources. 


'd 
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The  quadrilateral  size  is  crucial  in  obtaining  the  solution  of  the 
Fredholm  integral  equation  and  then  calculating  the  hydrodynamic  forces  and 
moments.  The  dimensions  are  dictated  by  the  propeller-induced  velocity 
which  is  of  oscillatory  character.  The  patch  dimensions  should  be  selected 
so  as  to  fit  at  least  four  quadrilaterals  into  one  cycle  of  spatial  oscil¬ 
lation  wherever  such  oscillations  are  present.  This  requirement,  however, 
increases  the  already  large  number  of  quadrilaterals  considerably  and, 
hence,  the  number  of  field  points  at  which  the  three  velocity  components  are 
to  be  evaluated  with  all  the  required  combinations  of  the  frequencies  of  the 
loading  and  propagation  functions. 

A  new  study  has  been  undertaken  at  Davidson  Laboratory  to  replace  the 
present  method,  which  represents  the  hull  by  a  source  distribution  and  re¬ 
quires  calculation  of  three  velocity  components  at  each  field  point  (null 
point),  by  a  procedure  which  represents  the  hull  by  a  doublet  distribution 
and  requires  the  calculation  of  only  one  function,  the  propeller  velocity 
potential.  This  method  permits  decreasing  the  mesh  size  of  the  quadrilateral: 
and  increasing  the  number  while  still  reducing  computer  time  and  is  shown  to 
be  a  considerable  improvement. 

E.  POINT  PRESSURES  ON  THE  HULL 

1 )  Due  to  the  Diffraction  Potential 

Having  obtained  a  ,  the  strength  of  the  sources  through  Eq.f1*),  the 
diffraction  potential  cp^  (designated  as  in  Eq.(3))  can  be  calculated, 
and  by  means  of  the  extended  Lagally  Theorem,  the  unsteady  hull  force  can 
be  determined  (see  Eq.(2i+),  Ref. 3).  Most  importantly,  the  pressure  field 
due  to  the  diffraction  potential  is  given  as 

/dtPD  StPD\ 

P0(X,Y,Z;t)  =  -Pf(sr  +  u  w~)  (24) 

If  it  is  assumed  that  cp^  has  simple  harmonic  time  dependence,  in  fact, 
considering  the  blade  frequency  component  only, 


iNClt 

e 


then  the  blade  frequency  pressure  due  to  the  diffraction  potential  is 

PD  -  -ipNfi  |9d|-  PU  u°  (25) 

where  uP  ,  the  corresponding  axial  component  of  the  velocity  is  determined 

A 
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as  a  by-product  of  the  process  of  inverting  the  integral  Eq.(4)  by  the 
Hess-Smith  numerical  code.  It  is  seen  that  the  convective  term  of  the 
blade  frequency  pressure  due  to  the  diffraction  potential  can  easily  be 
determined  from  the  available  information.  The  time-dependent  term  re¬ 
quires  an  explicit  expression  of  cp^. 

No  attempt  has  been  made  to  develop  a  computer  program  for  the  eval¬ 
uation  of  the  diffraction  potential  cp^  because  of  limited  funding.  The 
diffraction  potential  could  be  determined  through  an  integration  of  the 
vertical  components  of  the  velocity,  ,  along  a  vertical  line  from  a 
given  point  on  the  hull,  where  the  pressure  is  to  be  calculated,  to  the 
waterplane  (free  surface)  where  the  potential  is  zero  because  of  the  invoca¬ 
tion  of  the  high  frequency  Froude  number  condition. 


Through  this  vertical  integration 


S  n 

J  w"(X,Y,Z)dZ 


®S  dtp 

J  dz  =  cpD(ts)  -  cp°(CH)  -  -  9D(Ch) 


since  9°(CS)  =0  on  the  waterplane.  Thus,  the  diffraction  potential  cp  is 
given  by 

=  -  J  w°(x,Y,z)dz  '  (26) 

The  vertical  diffraction  velocities  required  in  this  integration  are  given 
by  successive  application  of  the  Hess-Smith  program  for  off-body  points^ 
which  may  be  particularized  for  points  located  inside  the  body.  Care 
should  be  taken  for  points  close  to  the  surface  of  the  body  where  there 
exists  a  sharp  gradient  of  the  velocity.  Another  difficulty  is  the  fact 
that  the  smooth  body  is  approximated  by  planar  quadrilaterals  so  that  the 
exact  location  of  the  body  surface  is  no  longer  precisely  known.  To  over - 
come  this  difficulty,  a  fine  mesh  size  should  be  used.  By  plotting  the 
results  of  the  calculated  velocities,  the  null  points  on  the  approximated 
body  will  clearly  be  indicated  by  the  presence  of  a  sharp  discontinuity  in 
velocity  at  the  points  across  the  boundary.  Thus,  the  ordinates  of  the  null 
points  and  the  range  of  integration  will  be  established. 

All  quantities  have  been  computed  in  non-dimensional  form,  using  the 
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model  speed  U  and  propeller  radius  rQ  (for  non-d imens ional i z i ng  the 
linear  dimensions).  The  non-dimensional  diffraction  velocity  potential 
(Eq.26)  is  g  i  ven  by 


D 


-Ur. 


(27) 


and  the  corresponding  expression  for  the  blade  frequency  pressure  (Eq.25) 
is  given  by 


+ipN2nn  rQ 


I 

C' 


W^dZ  1 


(28) 


with  the  substitution  Q  =  2TTn,  n  =  rps. 


2)  Free  Space  Pressure 

These  results  must  then  be  combined  with  pressures  emanating  from  the 
propeller  in  the  absence  of  the  hull  surface.  These  pressures  are  calculated 
utilizing  the  pressure  field  program^'0^  in  conjunction  with  the  blade  pres¬ 
sure  program  of  Reference  7 .  For  completeness  of  the  present  study,  the 
expressions  for  these  pressures  are  given  below: 

i  )  Free  space  pressure  due  to  propeller  loading 


P(*.r,cp;t)  =  J 


oo  n  max 

£  S  L.  -(p)e 
\=o  n=l 


i?vft  t  -imftt  ,,  /  % 

Z/  e  K  -  x.r.cpjp) 

m,  n 


m=-ce 


(29) 


where,  with  p,r,x,k  and  a  non-dimensional i zed  with  respect  to  propeller 
radius  rQ  ,  the  influence  or  spatial  propagation  function  is 


,  \  Ni  im(a'-cp) 

V^lx.r^p)  -+  — « 

o 

•I*  J  klm(kp)Kn,(kr)teik(x-CT,/a)A("> 

o 


+  V  j  lm(kp)K„(kr)[.lk<«-0,/a)A(5) 
P  o 


((m-|)6b)-e- 
((m“  l)eb)+e’ 


k(x-oVa)A(R)(m+^)6b)]dk 
k(x-o'/a)/n) (m+  ^)6b)]dk| 


(30) 
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where  for  m^O  and  k=0  ,  the  integrand  is  equal  to 


i  1 m I  ,-s 

A(n)(m0b)  for  P  <  r 


and 


i  -  lmf  f-\ 

- —  (7)  A  n  (m6b)  for  r  <  p 


p2  Iml  vp/ 

When  m=0  ,  Eq.(30)  becomes 

Ko,R(x»r*Cp;p)  =  +  ^7;  i23'  I  k,o^kp^Ko^kr^ 

o  0 

•  imag.  P.rt[.lk<,-<’,/*)A(">(-  J  8b>]  } 

In  this  case  when  k=0  ,  the  integrand  is  zero. 


(30 


In  computing  the  k-integral  of  (30)  or  (31)  numerically,  it  must  be 
realized  that  the  integration  interval  Ak  depends  mainly  on  the  range  of 
x.  The  product  of  the  modified  Bessel  functions  decreases  monoton ical ly 

with  increasing  k,  and  the  A  function,  although  oscillatory,  is  damped  out 

with  increasing  k.  However,  exp(±ikx)  is  undamped  oscillatory  and  care 
must  be  taken  to  have  sufficient  ordinates  in  a  cycle. 

Using  the  following  designation  for  simplicity 
8  p.+0 

L\  n  km  n  =  S  l\  n(pi)  f  Km  n  (* .  r , <pi  p)  d p 

A ,  n  m  j  n  •  ■  A  j  n  i  j  fn  j  n 

1=1  P..p 


where  20  is  the  spanwise  strip  length,  then  the  pressure  can  be  expressed  as 


n  max  «  ® 


P(xfr,<p;t)  «  E  E  E  X,  -  k  - 


i  ( X-m)  fit 


n=l  X=o  m=-a 


X ,  n  m ,  n 


(32) 


where  only  those  terms  must  be  summed  for  which  X-m=iN  and  X-m“-£N  , 
X  *  0, 1 ,2,  ... 

In  the  steady-state  case  (mean  value)  when  £=0  ,  then 
n  max 


Pc(x,  r  ,cp)  «  E  E  l  -  k  - 

n-1  X-o  *’n  K,n 


(33) 


22 


R-2140 


In  the  unsteady  case  (£  >  0)  the  pressure  given  by  Eq .(32)  can  be  written  as 


P(x,r,cp;t)  =  R.P. 


2 

£=o 


n  max 
2 
n=l 


=  [ 

X=o 


L  - 

X,  n 


kX-XN,n  e 


-UN6 


A,n 


kX+£N,ne 


+  i  £N0* 


=  R.P. 


®  n  max 
2  _2 
£=o  n=l 


sK.n 

X=o 


CX-£N, n  +conJ*  *L^\,n  kX+£N,n]  }e 

(34) 


UN6 


where  Qt  has  been  replaced  by  -0,  the  blade  angular  position  (see  Figure  3). 

The  complex  amplitude  of  the  £th  blade  frequency  harmonic  pressure 
i  s 


n  max  »  .  -  ... 

P£N  =  P.  kX-£N,n  +  COnj’ ‘  lA,n  kX+£N,nJ  / 

n  —  I  a — o 

It  may  be  observed  that  the  pressure  field  component  at  frequency  £N 
generated  by  an  operating  propeller  in  nonuniform  inflow  is  made  up  of 
linear  combinations  of  terms  proportional  to  blade  loadings  at  all  significant 
wake  harmonics  X  propagated  by  spatial  functions  of  harmonic  orders  m=X-£N 
and  m  =  X+£N. 

It  is  interesting  to  note  that  the  propeller-induced  pressure  field 
in  the  uniform  inflow  case  is  remarkably  different  from  that  in  the  non- 
uniform  case.  In  the  uniform  case,  only  the  zero  harmonic  of  loading  exists, 
i.e.,  X=0,  and  hence  the  corresponding  pressure  field  is  given  by 


1)  For  the  mean  (mean  value,  £= 0,  therefore  m=0) 


n  max 


P  (x,r)  =  2  £  -  k  - 

ov  '  -  .  o,n  o,n 

n=l 


(36) 


2)  For  the  £Nth  harmonic 


P£N(x,r,cp) 


n  max 
2  2 
n=l 


^o,n 


k-£N,n 


(37) 


One  should  appreciate  that  these  amplitudes  are  strong  near  the  propeller,  but 
decay  extremely  rapidly  with  distance  from  the  propeller. 
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i  i )  Pressure  field  due  to  propeller  thickness 

The  effect  of  finite  blade  thickness  is  to  produce  a  symmetrical 
disturbance  of  the  fluid  on  each  side  of  the  surface  of  symmetry  through 
the  blade  sections.  The  thickness  distribution  of  the  blade  section  is 
approximated  by  a  source-sink  distribution  smeared  over  the  section  with 
strength  determined  by  the  "thin  body"  approximation. 

Reference  10  gives  the  pressure  due  to  blade  thickness  at  multiples 
of  blade  frequency  as 


<%>*<'<■» 


■  J  uF(u)[G(u,  i)  +  G(-u,  jfc)  jdudpj-  (38) 


where 


,  r  .  ueb  u6b  u0b-|  /« 
F(u)  =Lsin-i - TCOS— J/U 


G(u,£)  =  l^N(lu+a£Nlp)K^N(lu+a£NI  r)e'U^X  °  ^ 


for  p<  r  ,  and  the  negative  sign  of  L  is  taken  into  account  in  the  com¬ 
puter  program. 

The  contribution  to  the  mean  or  zero  harmonic  is  given  by 


<PA 


t  (  P)  08 

TfpT  |  uF(u)lo(up)Ko(ur) 


•  cos[u(x-a'/a)]dudpj- 


(39) 


A  numerical  procedure  for  the  evaluation  of  the  pressure  field  due  to 
propeller  loading  and  thickness  designated  as  PFDTLOG  is  coded  in  Reference  10, 
provided  the  spanwise  loading  along  the  blade  at  the  desired  frequencies  is 
known  through  the  use  of  PPEXACT^ 

Combining  Eqs.(28),  (34)  and  (38),  the  blade  frequency  pressure  at  any 
point  on  the  hull  surface  can  be  determined  due  to  the  direct  propeller 
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pressure  (i.e.,  propeller  fanning  effect)  and  due  to  the  diffraction  on 
account  of  the  presence  of  the  hull. 

F.  CORRELATIONS  OF  THEORY  WITH  EXPERIMENTS 

In  order  to  establish  the  reliability  of  the  procedure  for  predicting 
the  mean  and  blade-frequency  pressure,  the  theoretical  calculations  have 

(4) 

been  compared  with  an  extensive  set  of  measurements  made  in  Norway.  Huse 
measured  propeller-induced  pressures  at  34  points  on  the  after-body  of  a 
20-foot  tanker  model  (l;40  scale)  of  block  coefficient  0.80  in  the  Norwegian 
Ship  Model  Experiment  Tank.  The  main  particulars  of  the  ship  and  model  are 
given  in  Table  1,  and  of  the  propeller,  in  Table  2. 


TABLE  1 

PARTICULARS  OF  MODEL  699 

Prototype  Model 

Length  between  perpendiculars  250  m  20.5' 

Breadth  42  m  3.45' 

Draft  15.24m  1.25' 

Block  coefficient  0.80  0.80 

Ship  speed  17  kt  4.54  ft/sec 


TABLE  2 

CHARACTERISTICS  OF  PROPELLER  TR00ST  B-4  SERIES 


Number  of  blades 
EAR 

Diameter,  D,  ft 

P/D 

RPM 

Ship  speed,  U,  ft/sec 
Advance  ratio,  J  =  U/nD 


4 

0.55 

0.671 

0.56 

680 

4.54 

0.6 

5.236 
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Relevant  geometric  characteristics  of  the  Troost  B-4  propel ler, namely, 
the  ratios  of  maximum  camber  to  chord  length,  m^/c,  maximum  thickness  to 
chord,  tQ/c,  of  leading  edge  radius  to  chord,  pQ,  of  chord  to  diameter,  c/D, 
and  of  the  pitch  to  diameter,  p/D,  at  eight  radial  positions,  are  given  in 
Table  3. 


TABLE  3 

PERTINENT  GEOMETRIC  CHARACTERISTICS  OF  TROOST  B-4  PROPELLER 


Radi  us 

m  /c 

X 

t  /c 
o 

po 

c/D 

p/D 

0.25 

0.0275 

0.1700 

0.0200 

0.223 

0.56 

0.35 

0.0260 

0.1300 

0.0135 

0.256 

0.56 

0.45 

0.0260 

0.1030 

0.0092 

0.282 

0.56 

0.55 

0.0254 

0.0820 

0.0058 

0.301 

0.56 

0.65 

0.0235 

0.0620 

0.0029 

0.308 

0.56 

0.75 

0.0213 

0.0488 

0.0012 

0.298 

0.56 

0.85 

0.0187 

0.0391 

0.0005 

0.259 

0.56 

0.95 

0.0137 

0.0313 

0.0002 

0.171 

0.56 

The  Troost 

B-4  propel le 

r  operates  in  the  wake 

of  Model 

699.  The  har 

content  of 

the  non-dimensional 

normal  velocity  on  the 

blade,  at  8 

stant  radial  positions 

between 

prope 1 ler  hub 

and  tip, 

is  given  in 

Table  4. 


TABLE  4 

HARMONIC  CONTENT  OF  THE  NON-DIMENSIONAL  NORMAL  VELOCITIES 
TO  THE  PROPELLER  BLADE  AT  VARIOUS  RADIAL  POSITIONS 


v(q)/  m  (Ref.  12) 
VN  ' 


q  * 

_ 9_! 

=  1 

q  =■ 

2 

q  s 

•  3 

q  = 

4 

q  *= 

5 

r/ro 

D 

D 

n 

D 

R 

n 

B 

B 

B 

B 

R 

B 

.25 

-.508 

0 

BviiLI 

BB 

+.067 

0 

-.021 

-.011 

+  .015 

MM 

.35 

-.557 

0 

BBS 

■91 

191 1'>] 

■Hm 

Bd 

+.019 

+.012 

m 

-.036 

+.021 

RSpf 

B0 

.45 

-.580 

0 

-.043 

tti 

+.008 

-.028 

M 

-.056 

+.023 

-.004 

.55 

-.592 

0 

g||| 

-.018 

-.001 

+.012 

-.043 

-.005 

-.065 

+.022 

-.023 

+.010 

.65 

-.604 

0 

-.002 

0 

+  .014 

-.027 

+.002 

-.050 

+.017 

-.029 

+.013 

.75 

-.574 

0 

+.036 

+.008 

-.055 

+  .005 

-.009 

+.007 

-.029 

+.016 

-.034 

+.009 

.85 

-.509 

0 

+.020 

am 

-.116 

+  .017 

-.010 

+.007 

-.038 

+.016 

-.024 

+.004 

.95 

0 

0 

BS 

-.209 

+  .032 

+  .063 

-.008 

-.083 

+.016 

-.008 

+  .002 

With  the  geometric  propeller  information  and  harmonic  content  of  the 
inflow  field  in  which  the  propeller  operates,  the  steady  and  time-dependent 
blade  loading  and  corresponding  hydrodynamic  forces  and  moments  can  be  de¬ 
termined  by  utilizing  the  numerical  procedure^  developed  and  programmed 
for  CDC-6600  or  Cyber  176  high-speed  digital  computer. 

The  results  of  the  calculations  for  the  steady  and  unsteady  thrust, 
torque  and  bearing  forces  and  moments  are  given  in  Table  5  (see  Figure  4) 
together  with  the  results  for  steady  state  thrust  and  torque  obtained  from 
Troost-charts  for  comparison. 

TABLE  5 
SUMMARY  OF 

HYDRODYNAMIC  FORCES  AND  MOMENTS  FOR  TROOST  B-4  PROPELLER 
IN  WAKE  OF  MODEL  699 


Steady  State 

Unsteady  (Blade  Frequency) 

Calculations 

Calculations 

i<T  =  0.2278 

0.173 

K  =  0.01601 

Kq  =  0.0244 

0.0178 

Kq  =  0.001428 

Kc  =  0.00145 
hy 

Kf  =  0.00104 

ry 

\L  =  0.00140 
vy 

Kq  =  0.00125 

Krr  =-0.00174 
z 

Kf  =  0.00102 
rz 

Kqz  =-0.00048 

Kq2  =  0.00158 

For  the  purpose  of  obtaining  the  solution  of  integral  equation  (4), 
the  hull  surface  and  its  positive  reflection  in  the  free  surface,  are  sub¬ 
divided  into  392  quadrilaterals  as  shown  in  Figures  5  and  6.  The  coded 
numerical  procedure  developed  by  Hess-Smith  is  applied,  considering  the 
following  onset  flows;  a  uniform  steady  motion  of  the  hull  and  its  image, 
yielding  the  source  density  o  and  the  real  and  imaginary  parts  of  the 
propeller-induced  onset  flows  which  determine  the  corresponding  strengths 

a  and  a  ,  respectively. 

R  I 

That  the  size  of  the  quadrilateral  is  a  crucial  factor  not  only  in 
calculating  the  hydrodynamic  forces  and  moments  but  also  in  evaluating  the 
point  pressure  and  velocity  (both  oscillatory  in  character)  on  the  hull  has 
been  previously  discussed.  The  selection  of  proper  dimensions  for  the  various 
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patches  considerably  increases  the  number  of  the  required  quadrilaterals. 

This  coupled  with  the  complications  arising  in  the  calculation  of  the 
propeller-induced  velocity  components  makes  the  numerical  approach  expensive 
and  rather  impractical. 

The  approach  used  in  this  study  has  therefore  not  been  pursued  further. 

(13} 

As  mentioned  earlier,  a  new  study'  '  has  been  initiated  which  represents 
the  hull  by  a  doublet  distribution  whose  strength  is  calculated  by  the  same 
numerical  procedure  for  the  inversion  of  the  integral  equation,  knowing  the 
propeller  velocity  potential  (a  scalar  function).  This  method  permits  de¬ 
creasing  the  mesh-size  of  the  quadrilaterals  and  increasing  the  number 
while  still  reducing  the  computer  time,  and  has  shown  considerable  improve¬ 
ment. 

However,  for  this  report  the  calculations  have  been  performed  using 
the  quadrilateral  arrangement  shown  in  Figures  5  and  6.  Out  of  392  quadri¬ 
laterals,  142  located  at  the  forward  part  of  the  hull  and  away  from  the 
propeller  are  deactivated  since  the  propeller-induced  velocities  there  are 
negl ibi ly  smal 1 . 

The  procedure  previously  developed  in  section  E  is  applied  to  calcu¬ 
late  the  point  pressures  at  the  locations  of  the  transducers  shown  in  Table  6, 
which  are  referred  to  the  propeller  coordinate  system  (see  Figure  7).  The 
free  space  pressure  (see  E-2)  is  determined  by  means  of  the  numerical  pro¬ 
cedure  coded  in  Reference  10  as  PFOTLDG  and  the  pressure  to  the  diffraction 
potential  is  determined  according  to  the  procedure  discussed  in  section  E-l. 
The  results  of  each  are  presented  in  Table  7  as  are  the  total  pressures  with 
their  corresponding  magnitudes  and  phases  referred  to  the  hull  coordinate 
system. 

The  calculations  are  performed  for  a  propel ler-skeg  clearance  of 
b/D  *  0.235  (see  Figure  3.6  of  Reference  4).  It  should  be  noted  that  the 
transducer  coordinates  are  determined  from  the  values  of  §/D  of  Table  4.1  of 
Reference  4  as  follows: 

a)  For  points  on  the  transom  of  the  hull 

x/D  ■  (5/D  given  in  Table  4.1  at  b/D  *  0.48)  +  (0.48  -  0.235) 

b)  For  points  on  the  hull  skeg 

x/D  ■  (5/D  given  in  Table  4.1  at  b/D  ■  0.186)-  (0.235-  0.186) 
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The  compaiison  between  predicted  and  measured  magnitudes  and  phases 
is  presented  in  Table  8a,  and  graphically  exhibited  in  Figures  9,  10  and  11  at 
various  hull  cross-sections  (see  Figure  8).  The  cross-section  F-F  includes 
also  point  pressures  from  neighboring  sections,  a  fact  which  explains  the 
abrupt  change  in  the  pressure  distribution  at  some  points. 

The  results  in  Table  8  are  presented  in  pressure  coefficient  Cp  form 
by  defining 


C 


P 


P 

pn2D2 


where  p  =  magnitude  point  pressure  on  hull 
p  *  fluid  dens i ty 
D  =  propeller  diameter 
n  =  rps 

The  ratio  of  Cp  measured  to  Cp  calculated  values  indicates  the  degree  of 

correlation  between  experiment  and  theory.  Note  that  the  phase  angles  are 

referred  to  field  point  angular  position  so  that  the  phase  angle  cp  as  de- 

c 

fined  in  Reference  4  is  given  by 

cpc  =  270°  -  4e  +  cp^ 


where 


-1  I 


cp^  =  phase  angle  obtained  through  the  program  as  tan  '  —■ 
e  =  the  angular  position  of  the  field  point  given  in  Table  6. 


As  was  mentioned  earlier,  the  present  method  where  the  hull  is  represented 
by  a  distribution  of  sources  smeared  over  the  hull  surface  has  been  found  un¬ 
economical  and  has  been  replaced  by  a  new  method  where  the  hull  surface  is  repre¬ 
sented  by  a  doublet  distribution.  Calculations  were  performed  with  the  same 
number  of  quadrilaterals  giving  a  detailed  pressure  distribution  which  is 
graphically  exhibited  in  Figure  9,  10  and  11.  On  those  graphs,  both  the  ex¬ 
perimental  results  and  the  results  using  the  source  distribution  are  presented 
as  well. 

A  correlation  of  theoretical  results  obtained  through  the  "doublet- 
approach"  with  those  of  experiments  is  presented  in  Table  8b.  It  would  be 
difficult  to  argue  the  superiority  of  the  "doublet-approach"  over  the 
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"source-distribution"  when  the  comparison  is  limited  to  the  magnitudes 
of  the  point  pressures  only.  When  the  theoretical  phase  angles  are  com¬ 
pared  with  those  of  the  corresponding  experiments,  it  is  definitely  con¬ 
cluded  that  the  "doublet-approach"  is  superior.  This  resul t,  coupled  wi th 
the  fact  that  the  "doublet-distribution"  is  more  economical  by  yielding 
greater  accuracy  without  using  excessive  computing  time,  demonstrates  the 
advisability  of  the  "doublet-approach"  rather  than  the  "source-distribution." 

In  Tables  9a  and  9b,  the  values  of  the  so-called  "reflection  coef- 
ficents"  are  presented  as  results  of  theoretical  calculations,  S^.  and 
a  combination  of  measurements  and  calculations,  Sj.  Table  9a  presents  the 
results  of  calculations  obtained  by  means  of  the  "source-approach"  whereas 
Table  9b  presents  the  results  obtained  by  the  "doublet-approach." 

(M 

The  "reflection  coefficient"  defined  by  Huse'  as  the  "solid 
boundary  factor"  actually  represents  the  "image  effect"  due  to  the  pres¬ 
ence  of  the  rigid  boundary.  It  is  known  that  for  a  rigid  boundary  of  a 
plane  surface  of  an  infinite  extent,  the  reflection  coefficient  yields  a 
value  of  S=Z.O  .  For  a  hull  surface  with  a  combination  of  convex,  concave, 
and  skeg-type  surfaces,  it  is  difficult  to  assume  any  suitable  value  before¬ 
hand.  It  is  certainly  true  that  the  approximation  commonly  used  by  doubling 
the  free-space  pressure  is  valid  in  cases  with  a  flat,  or  nearly  flat  bottom 
stern  in  the  way  of  the  propeller;  however,  there  are  many  other  instances  in 
which  the  approximation  is  very  rough  and  therefore  unacceptable.  None  of 
the  results  of  Tables  9a  and  9b  shows  any  definite  trend  in  this  direction. 
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TABLE  6 

TRANSDUCER'S  LOCATION  IN  TERMS  OF  PROPELLER  COORDINATES 
WITH  PROPELLER-SKEG  CLEARANCE  AT  b/D  =  0.235 


T  ransducer 

x/r0 

r/ro 

cp,  radians 

or  e,  deg 

C1 

T  ransom 

0.194 

1.514 

0 

0 

'2 

T  ransom 

-0.168 

1.420 

0 

0 

Pi 

Transom 

1.758 

2.040 

0.197 

11.3 

p  2 

T  ransom 

0.976 

1.860 

0.216 

12.4 

p3 

T  ransom 

0.194 

1 .682 

0.241 

13.8 

p4 

T  ransom 

-0.168 

1.600 

0.253 

14.5 

p5 

T  ransom 

-0.538 

1.464 

0.278 

15.9 

p7 

T  ransom 

0.194 

1.970 

0.396 

22.7 

p8 

T  ransom 

0.194 

2.300 

0.511 

29.3 

Pq 

T  ransom 

0.194 

2.652 

0.588 

33.7 

P10 

Skeg 

-0.684  (b) 

1.020 

0.153 

8.8 

pn 

Skeg 

-0.722 

0.728 

0.202 

11.6 

P12 

Skeg 

-1.016 

0.772 

0.389 

22.3 

P13 

Skeg 

-0.636 

0.454 

0.328 

18.8 

pl4 

Skeg 

-0.576 

0.446 

2.874 

164.7 

P15 

Skeg 

-0.748 

0.760 

3.00 

171.5 

Pl6 

Skeg 

-1.046 

0.774 

2.899 

166.1 

S1 

Transom 

1.758 

2.040 

-0.197 

-  11.3 

=2 

T  ransom 

0.976 

1.860 

-0.216 

-  12.4 

S3 

T  ransom 

0.194 

1.682 

-0.241 

-  13.8 

54 

T  ransom 

-0.168 

1.600 

-0.253 

-  14.5 

55 

Transom 

-0.538 

1 .464 

-0.278 

-  15.9 

S7 

T  ransom 

0.194 

1.970 

-0.396 

-  22.7 

s8 

Transom 

0.194 

2.300 

-0.511 

-  29.3 

5,0 

Skeg 

-0.684 

1.020 

-0.153 

-  8.8 

S11 

Skeg 

-0.722 

0.728 

-0.202 

-  11.6 

S12 

Skeg 

-1.016 

0.772 

-0.389 

-  22.3 

S13 

Skeg 

-0.636 

0.454 

-0.328 

-  18.8 

5,4 

Skeg 

-0.576 

0.446 

-2.874 

-164.7 

S 15 

Skeg 

-0.748 

0.760 

-3.00 

-171.9 

Sl6 

Skeg 

-1.046 

0.774 

-2.899 

-166.1 
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TABLE  7 

POINT  PRESSURES  AT  THE  TRANSDUCER'S  LOCATION 
DUE  TO  PROPELLER  AND  DUE  TO  DIFFRACTION  POTENTIAL  (Presence  of  the  Hull) 


T  ransd 

Quad. 

; 

THEORETICAL  POINT 

PRESSURE 

Due  to  Propeller 

Due  to  Diffraction 

Total  Pressure  I 

Real 

Imaginary 

Real  Imaginary 

Real 

Imaginary 

Mag 

9 

C1  > 

k  12,13 

0.26351 

0.21697 

-0.38404  0.20936 

-0.12053 

0.42633 

0.443 

106 

‘2 

14 

-0.39576 

0.28795 

-1.07607  0.61195 

-1.47183 

0.89990 

1.725 

149 

P1 

17 

0.06818 

0.05127 

0.01346  0.02186 

0.08164 

0.07313 

0.110 

42 

0.07047 

0.04393 

0.02139  0.03280 

0.09186 

0.07673 

0.120 

40 

P, 

19 

0.10533 

0.09185 

-0.00192  0.01741 

0.10341 

0.10926 

0.150 

47 

^2 

0.11548 

0.02331 

0.03533  0.07155 

0.15081 

0.09486 

0.178 

32 

Pr 

22 

0.00640 

0.15652 

-0.21569  -0.05289 

-0.20929 

0.10363 

0.234 

154 

s3  O 

0.18377 

0.01966 

-0.02008  0.27974 

0.16369 

0.29940 

0.341 

61 

D  CO 

p4  5 

25 

-0.25249 

-0.03248 

-0.66811  -0.37877 

-0.92060 

-0.41125 

1.008 

-156 

<*  2 

0.03646 

0.16795 

-0.15068  0.72955 

-0.11422 

0.89750 

0.905 

97 

p5 

28 

-0.35090 

-0.22753 

-0.52340  -0.62954 

-0.87430 

-0.85707 

1 .224 

-136 

S5 

-0.12372 

0.17732 

0.10715  0.75310 

-0.01657 

0.93042 

0.931 

91 

p6 

31 

0.05339 

0.06329 

-0.00346  -0.01182 

0.04993 

0.05147 

0.070 

46 

p7 

34 

-0.02829 

0.04017 

-0.06663  -0.14383 

-0.09492 

-0.10366 

0.141 

-132 

S7 

0.07687 

-0.01068 

0.11031  0.14798 

O.I8718 

0.13730 

0.232 

36 

p8 

33 

-0.01605 

0.00246 

-0.02182  -0.09106 

-0.03787 

-0.08860 

0.096 

-113 

s8  , 

0.03431 

-0.00916 

0.07820  0.06949 

0.11251 

0.06033 

0.128 

28 

p9 

32 

-O.OO676 

-0.00767 

-0.00564  -0.04532 

-0.01240 

-0.05299 

0.054 

-103 

P,0  /| 

fs, 

44 

-0.85609 

-0.28651 

-0.67435  -0.49409 

-1.53044 

-0.78060 

1.718 

-153 

S10 

-0.55158 

0.31289 

0.0^075  0.75641 

-0.48083 

1.06930 

1.172 

114 

P11 

45 

-0.56181 

-0.37761 

-0.40035  -0.63409 

-0.96216 

-1.01170 

1.396 

-133 

S11 

-0.25675 

0.12826 

0.07047  0.71052 

-0.18628 

0.83878 

0.859 

103 

P1 2 

58 

-0.10162 

-0.43533 

-0.08001  -0.02526 

-0.18163 

-0.46059 

0.495 

-111 

S12 

-0.33155 

0.06926 

0.14401  0.24039 

-0.18754 

0.30965 

0.362 

121 

P13 

46 

-0.48863 

-0.75591 

-0.03826  -0.1/861 

-0.52689 

-0.93452 

1.073 

-119 

sn 

D 

-0.42611 

-0.02570 

-0.12649  0.27476 

-0.55260 

0.24906 

0.606 

156 

p 

P14  ^ 

48 

-0.65826 

0.45380 

-0.90732  0.07214 

-1.56558 

0.52594 

1.652 

161 

SI4 

-0.71043 

-0.72396 

-0.47814  0.01274 

-1.18857 

-0.71122 

1.385 

-149 

P15 

53 

-0.51649 

0.35695 

-0.18960  0.13821 

-0.70609 

0.49516 

0.862 

145 

S15 

-0.74848 

-0.10680 

-0.15475  -0.25284 

-0.90323 

-0.35964 

0.972 

-158 

P16 

62 

-0.39088 

-0.05751 

0.01550  0.06337 

-0.37538 

0.00586 

0.375 

179 

S16 

-0.11371 

-0.26772 

-0.11225  -0.30656 

-0.22596 

-0.57428 

0.617 

-111 

P17 

63 

-0.26621 

0.15803 

s,7s 

* 

-0.36241 

0.05653 
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TABLE  8a 

CORRELATION  OF  THEORETICAL  RESULTS  WITH  MEASUREMENTS 
SOURCE-APPROACH 


1 

COMPARISON  OF  THEORETICAL 

AND  EXPERIMENTAL  RESULTS  "1 

Theoretical 

Experimental 

KM 

ex  th 
9  -9 

c  th 

UP 

'Pc.  deg 

Upper 

Boundary 

Lower 
Bounda  ry 

c  ex 

LP 

9c>de9 

c  th 

LP 

deg 

0.00396 

16 

0.00498 

0.00456 

0.0048 

236 

1.21 

220 

0.015^1 

59 

0.0157 

0.0155 

0.0156 

268 

1.01 

209 

0.00098 

- 

0.00076 

- 

0.78 

0.00107 

- 

0.00106 

0.00054 

0.0009 

- 

0.84 

0.00134 

- 

0.00084 

0.00036 

0.00060 

- 

0.45 

0.00159 

- 

0.00275 

0.0021 

0.0024 

- 

1.50 

0.00209 

8 

0.0043 

0.0041 

0.0042 

263 

2.00 

'255 

0.00305 

27 

0.0063 

0.00610 

0.0062 

287 

2.03 

260 

0.00901 

56 

0.0115 

0.0111 

0.0113 

275 

1.25 

219 

0.00808 

65 

0.0103 

0.0099 

0.0101 

297 

1.25 

232 

0.01094 

71 

0.0105 

0.0095- 

0.0100 

277 

0.91 

206 

0.00831 

65 

0.0086 

0.0082 

0.0084 

298 

1.01 

233 

0.00064 

- 

- 

- 

- 

0.00126 

47 

0.0023 

0.0021 

0.0022 

250 

1.75 

203 

0.00207 

37 

0.0030 

0.0027 

0.0028 

264 

1.35 

227 

0.00086 

40 

0.0012 

0.0008 

0.0010 

245 

1.16 

205 

0.001 14 

55 

0.00155 

0.00125 

0.0014 

310 

1.23 

255 

0.00049 

- 

- 

- 

0.01535 

82 

0.0140 

0.0130 

0.0135 

298 

0 .88 

216 

0.01047 

59 

0.01630 

0.01530 

0.0156 

268 

1.49 

209 

0.01247 

90 

0.01300 

0.01200 

0.0125 

291 

1 .002 

201 

0.00767 

59 

0.01410 

0.0131 

0.0135 

CM 

1.76 

214 

0.00442 

69 

0.00470 

0.00430 

0.0045 

273 

1.02 

204 

0.00323 

120 

0.00520 

0.0050 

0.0051 

322 

1.58 

202 

0.00958 

75 

0.0078 

0.0068 

0.0075 

282 

0.78 

207 

0.00541 

141 

0.0055 

0.0045 

0.0050 

282 

0.93 

141 

0.01475 

133 

0.0100 

0.0095 

0.0100 

306 

0.68 

173 

0.01237 

60 

0.0095 

0.0085 

0.0090 

278 

0.73 

218 

0.00770 

87 

0.0100 

0.0090 

0.0095 

273 

1.23 

186 

0.00868 

79 

0.0090 

0.0080 

0.0085 

CO 

CO 

CM 

0.98 

204 

0.00335 

145 

0.0039 

0.0031 

0.0035 

288 

1.04 

143 

0.00551 

103 

0.0042 

0.0032 

0.0037 

256 

0.67 

153 
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TABLE  8b 

CORRELATION  OF  THEORETICAL  RESULTS  WITH  MEASUREMENTS 
DOUBLET -APPROACH 


Theoretical 

Experimental 

Cpex 

r  th 

LP 

<pc,deg 

Upper 

Boundary 

Lower 

Boundary 

c  ex 

LP 

9c » de9 

c  tb 

LP 

0.00559 

228 

0.00498 

0.00456 

0.0048 

236 

0.86 

0,02191 

261 

0.0157 

0.0155 

0.0156 

268 

0.71 

0.00105 

244 

- 

- 

0.00076 

- 

0.72 

0.00089 

291 

0.00106 

0.00054 

0.0009 

_.L 

1.00 

0.00032 

240 

0.00084 

0.00036 

0.00060 

- 

1.87 

0.00073 

308 

0.00275 

0.0021 

0.0024 

- 

3.28 

0.00499 

230 

0.0043 

0.0041 

0.0042 

263 

0.84 

0.00464 

239 

0.0063 

0.0061 

0.0062 

287 

1.34 

0.01445 

240 

0.0115 

0.0111 

0.0113 

275 

0.78 

0.01190 

256 

0.0103 

0.0099 

0.0101 

297 

0.86 

0.01741 

248 

0.0105 

0.0095 

0.0100 

277 

0.57 

0.01330 

269 

0.0086 

0.0082 

0.0084 

298 

0.63 

0.00179 

- 

- 

- 

0.00304 

226 

0.0023 

0.0021 

0.0022 

250 

0.72 

0.00156 

260 

0.0030 

0.0027 

0,0028 

264 

1.79 

0.00230 

199 

0.0012 

0.0008 

0.0010 

245 

0.43 

0.00047 

340 

0.00155 

0.00125 

0.0014 

310 

2.98 

0.00267 

173 

- 

- 

- 

- 

0.01903 

294 

0.0140 

0.0130 

0.0135 

298 

0.71 

0.01330 

251 

0.01630 

0.01530 

0.0156 

268 

1.17 

0.01950 

284 

0.01300 

0.01200 

0.0125 

291 

0.64 

0.01238 

271 

0.01410 

0.0131 

0.0135 

273 

■ 

tr 

O' 

00 

0 

0 

c 

241 

0.00470 

0.0043 

0.0045 

273 

0.00262 

309 

0.0052 

0.0050 

0.0051 

322 

1.95 

0.01202 

264 

0.0078 

0.0068 

0.0075 

282 

0.62 

0.00642 

350 

0.0055 

0.0045 

0.0050 

282 

0.78 

0.01261 

314 

0.0100 

0.0095 

0.0100 

306 

0.79 

0.01490 

263 

0.0095 

0.0085 

0.0090 

278 

0.60 

0.01 180 

274 

0.0100 

0.0090 

0.0095 

273 

0.81 

0.01390 

m 

0.0090 

0.0080 

0.0085 

283 

0.61 

0.00434 

262 

0.0039 

0.0031 

0.0035 

288 

0.81 

0.00603 

254 

0.0042 

0.0032 

0.0037 

256 

0.61 

ex  th 
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TABLE  9a 

REFLECTION  COEFFICIENTS  (Solid  Boundary  Factor) 
SOURCE-APPROACH 


Cp  Experiment  Cp  Total 

C  =  — -  .I  ■  —  -  — —  —  £  s  .  ■  -  -  ■  ■  ■  I 

1  Cp  Free-Space  Pressure  ’  2  Cp  Free-Space  Pressurel Theory 


T  rans 

Quad. 

THEORY 

EXPER. 

CP 

r  exp 

sr  fn, 

LP 

CptK 

s2=  fsp 

r 

FREE-SPACE  PRESSURE 

Due  to  Propeller 

Real  Imaginary  Magni.  Cp 

Total 

Pressure 

CP 

q  / 

k  12,13 

0.26351  0.21697 

0.34134  0.0030 

0.00396 

0.0048 

1.60 

c2 

1*4 

-0.39576  0.28795 

0.48943  0.0044 

0.01541 

0.0156 

3.54 

Pi 

17 

0.06818  0.05127 

0.08531  0.00076 

0.00098 

0.000# 

1.00 

1.28 

q 

0.07047  0.04393 

0.08304  0.00074 

0.00107 

0.0009 

1.21 

1.44 

p2 

19 

0.10533  0.09185 

0.13975  0.00124 

0.00134 

0.0006 

0.50 

1.08 

S2  3 

- 

0.11548  0.02331 

0.11781  0.00105 

0.00159 

0.0024 

2.28 

1.51 

p3  S  22 

0.00640  0.15652 

0.15665  0.00140 

0.00209 

0.0042 

3.00 

1.49 

S3  | 

0.18377  0.01966 

0.18482  0.00165 

0.00305 

0.0062 

3.75 

1.84 

P4  ^  25 

-0.25249  -0.03248 

0.25457  0.00227 

0.00901 

0.01  13 

4.97 

3.96 

S4 

0.03646  0.16795 

0.17186  0.00153 

0.00808 

0.0 10] 

6.60 

5.28 

Ps 

28 

-0.35090  -0.22753 

0.41826  0.00374 

0.01094 

0.0J00 

2.67 

2.92 

S5 

-0.12372  0.17732 

0.21622  0.00193 

0.00831 

0.0084 

4.35 

4.30 

Pe 

31 

0.05339  0.06329 

0.08280  0.00074 

0.00064 

- 

- 

- 

Pv 

34 

-0.02829  0.04017 

0.04913  0.00044 

0.00126 

0.0022 

5.00 

2.86 

S7 

0.07687  -0.01068 

0.07761  0.00069 

0.00207 

0.0028 

4.05 

3.00 

Pe 

33 

-0.01605  0.00246 

0.01624  0.00014 

0.00086 

0.0010 

7.14 

6.14 

Se  > 

f 

0.03431  -0.00916 

0.03551  0.00032 

0.001 14 

0.0014 

4.37 

3.56 

Pe  1 

'  32 

-0.00676  -0.00767 

0.01022  0.00009 

0.00049 

•* 

5.44 

*4*4 

-0.85609  -0.28651 

0.90276  0.00806 

0.01535 

m 

1.67 

1 .90 

sio; 

-0.55158  0.31289 

0.63415  0.00566 

0.01047 

MfaEa 

2.75 

1.84 

Pit 

*45 

-0.56181  -0.37761 

0.67692  0.00604 

0.01247 

0.0125 

2.06 

2.06 

1 

-0.25675  0.12826 

0.28700  0.00256 

0.00767 

0.0135 

5.27 

2.99 

Pi  2 

58 

-0.10162  -0.43533 

0.44703  0.00399 

0.00442 

0.0045 

1.13 

1.11 

Sj  2 

-0.33155  0.06926 

0.33870  0.00302 

0.00323 

0.0051 

1.69 

1.06 

P,a 

46 

-0.48863  -0.75591 

0.90009  0.00804 

0.00958 

0.0075 

0.94 

1.19 

q  3  ts 

-0.42611  -0.02570 

0.42688  0.00381 

0.00541 

0.0050 

1.31 

1.42 

P14S  ^8 

-0.65826  0.45380 

0.94186  0.00841 

0.01475 

0.0100 

1.19 

1.75 

s*4  1/1 

-0.71043  -0.72396 

1.01431  0.00906 

0.01237 

0.0090 

0.99 

1.36 

PjS 

53 

-0.51649  0.35695 

0.62783  0.00561 

0.00770 

0.0095 

1.69 

1.14 

Sl5 

-0.74848  -0.10680 

0.75606  0.00675 

0.00868 

0.0085 

1.26 

1 .28 

P1  6 

62 

-0.39088  -0.05751 

0.39509  0.00353 

0.00335 

0.0035 

0.99 

0.95 

Sl6 

-0.11371  -0.26772 

0.29086  0.00259 

0.00551 

O.OO37 

1.43 

2.12 

P1  7 

.  63 

-0.26621  0.15803 

0.30958  0.00276 

- 

- 

- 

Sl7^ 

-0.36241  0.05653 

0.36679  0.00328 
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TABLE  3b 

REFLECTION  COEFFICIENTS  (Solid  Boundary  Factor) 
DOUBLET-APPROACH 


Cp  Experiment  Cp  Total 


'1  Cp  Free-Space 

Pressure  ’  b2 

Cp  Free-Space 

PressureiTheory 

1  THEORY 

EXPER. 

FREE-SPACE  PRESSURE 

Total 

rCXp 

Cpth 

Trans  Quad. 

Due  to  Prope 1 1 er 

Pressure 

CP 

s,=  — f 

V  “fSD 

Real  Imaginary  Magni.  Cp 

CP 

r 

1  r  TsP 
LP 

z  -  frsp 
cp 

Ci  > 

f  12,13 

0.26351  0.21697 

0.34134  0.0030 

0.0056 

0.0048 

1.60 

1.86 

Ca 

14 

-0.39576  0.28795 

0.48943  0.0044 

0.0219 

0.0156 

3.54 

4.98 

Pi 

17 

0.06818  0.05127 

O.O8531  O.OOO76 

0.0105 

O.OOOyfc 

1.00 

1.38 

Si 

0.07047  0.04393 

0.08304  0.00074 

0.00089 

0.0009 

1.21 

1 .20 

Pa 

19 

0.10533  0.09185 

0.13975.  0.00124 

0.00032 

0.0006 

0.48 

0.25 

Sa 

E 

0.11548  0.02331 

0.11781  0.00105 

0.00073 

0.0024 

2.28 

0.70 

Pa  S  22 

0.00640  0.15652 

0.15665  0.00140 

0.00499 

0.0042 

3.00 

3.56 

s3  i 

0.18377  0.01966 

0.18482  0.00165 

0.00464 

0.0062 

3.75 

2.81 

p*  i 

-  25 

-0.25249  -0.03248 

0.25457  0.00227 

0.01445 

0.01 13 

4.98 

6.36 

s4 

0.03646  0.16795 

0.17186  0.00153 

0.01180 

0.0101 

6.60 

7.71 

28 

-0.35090  -0.22753 

0.41826  0.00374 

0.01741 

0.0100 

2.67 

4.65 

Ss 

-0.12372  0.17732 

0.21622  0.00193 

0.01330 

0.0084 

4.35 

6.89 

Pe 

31 

0.05339  0.06329 

0.08280  0.00074 

0.00179 

- 

• 

P7 

34 

-0.02829  0.04017 

0.04913  0.00044 

0.0030*4 

0.0022 

5.00 

6.91 

S7 

0.07687  -0.01068 

0.07761  0.00069 

0.00156 

0.0028 

4.06 

2.26 

Pe 

33 

-0.01605  0.00246 

0.01624  0.00014 

0.00230 

0.0010 

7.14 

16.43 

Se  > 

r 

0.03431  -0.00916 

0.03551  0.00032 

0.00047 

0.0014 

4.37 

1.47 

Pa 

1  32 

-O.OO676  -O.OO767 

0.01022  0.00009 

0.00263 

- 

- 

29.67 

Pio  , 

44 

-0.85609  -0.28651 

0.90276  0.00806 

0.0190 

0.0135 

1.67 

2.36 

Sio" 

-0.55158  0.31289 

0.63415  0.00566 

0.0133 

0.0156 

2.76 

2.35 

P11 

45 

-0.56181  -0.37761 

0.67692  0.00604 

0.0195 

0.0125 

2.06 

3.23 

S11 

-0.25675  0.12826 

0.28700  0.00256 

0.0124 

0.0135 

5.27 

4.84 

Pi  2 

58 

-0.10162  -0.43533 

0.44703  0.00399 

0.0090 

0.0045 

1.13 

2.25 

Si  a 

-0.33155  0.06926 

0.33870  0.00302 

0.0026 

0.0051 

1.69 

0.86 

Pi  a 

46 

-0.48863  -0.75591 

0.90009  0.00804 

0.0120 

0.0075 

0.93 

1.49 

Si  3  0 

-0.42611  -0.02570 

0.42688  0.00381 

0.0064 

0.0050 

1.31 

1.68 

PuS  48 

-0.65826  0.45380 

0.94186  0.00841 

0.0126 

0.0100 

1.19 

1.50 

Si  4." 

> 

-0.71043  -0.72396 

1.01431  0.00906 

0.0149 

0.0090 

0.99 

1.64 

Pis 

53 

-0.51649  0.35695 

0.62783  0.00561 

0.0118 

0.0095 

1.69 

2.10 

Sis 

-0.74848  -0.10680 

0.75606  0.00675 

0.0139 

0.0085 

1.26 

2.06 

Pis 

62 

-0.39088  -0.05751 

0.39509  0.00353 

0.0043 

0.0035 

0.99 

1.22 

Si6 

-0.11371  -0.26772 

0.29086  0.00259 

0.0060 

0.0037 

1.43  ! 

2.31 

P17 

^  63 

-0.26621  0.15803 

0.30958  0.00276 

0.0038 

« 

_ 

1.38 

Si 7 

-0.36241  0.05653 

0.36679  0.00328 

0.0055 

• 

- 

1.68 
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G.  CONCLUSIONS 

The  previously  developed  method  for  calculating  the  propeller-induced 
ship  vibration^  has  been  extended  to  establish  a  numerical  approach  for 
evaluating  the  pressure  distribution  on  a  hull  surface  in  the  presence  of  the 
propeller  and  the  free  water  surface.  The  hull  is  represented  by  a  source 
distribution  smeared  over  the  hull  surface,  whose  strength  is  determined 
through  the  solution  of  a  Fredholm  integral  equation  of  the  2nd  kind  in  terms 
of  specified  propeller-induced  velocities.  To  carry  out  the  numerical  solu¬ 
tion  of  the  integral  equation,  the  hull  surface  is  subdivided  into  plane 
quadrilaterals  whose  sizes  and  numbers  are  very  crucial  in  calculating  the 
propeller-induced  hydrodynamic  forces  and  point  pressure  on  the  hull.  Selec¬ 
tion  of  the  proper  size  for  the  source  panels  increases  the  number  of  quadri¬ 
laterals  and  this,  coupled  with  complications  arising  from  calculations  of 
the  propeller-induced  velocity,  makes  this  numerical  approach  expensive  and 
rather  impractical.  In  fact,  the  need  for  accuracy  in  the  numerical  procedure 
is  more  pronounced  in  calculating  point  pressure  than  in  determining  the 
propeller-induced  vibratory  force  since  the  latter  represents  the  integrated 
effect  of  the  point  pressure  over  the  hull  surface.  Because  of  this,  a  new 
approach  has  been  developed  in  which  the  hull  is  represented  by  a  doublet 
distribution.  The  computing  time  is  reduced  to  almost  half  that  incurred  by 
the  present  method  although  the  mesh  size  has  been  decreased  and  the  number 
of  quadrilaterals  increased.  The  blade-frequency  point  pressure  is  determined 
by  superposing  the  free-space  pressure  emanating  from  the  propeller  and  the 
pressure  due  to  the  diffraction  potential.  The  former  is  calculated  accord¬ 
ing  to  the  theory  and  program  of  Reference  10  whereas  the  latter  is  determined 
through  an  approximate  method  developed  in  Section  E  of  the  present  study. 

A  series  of  calculations  was  performed  to  determine  point  pressures 
on  an  M-699  model  tanker  with  CD  =  0.80  for  which  extensive  measurements 

D 

were  made  at  the  Norwegian  Experimental  Tank.  The  correlation  between  experi¬ 
mental  values  and  those  predicted  by  this  proposed  method  shows  a  tolerable 
discrepancy  in  magnitude  of  pressure  signals  but  the  values  of  the  correspond¬ 
ing  phases  show  very  poor  agreement.  This  fact  is  attributed  mainly  to  the 
size  and  number  of  quadrilaterals  which  affect  the  solution  of  the  integral 
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equation  and,  consequently,  the  values  of  the  vibratory  forces  and  point  pres¬ 
sures  on  the  hull  surface,  all  of  which  have  not  yet  converged  to  their  final 
values.  This  approach  has  been  suspended  as  being  expensive  and  impractical.  £ 

The  new  method  suggested  in  which  the  hull  is  represented  by  a  doublet 
distribution  rather  than  a  source  distribution  is  a  decided  improvement.  In 
fact,  the  doublet  distribution  is  superior  to  the  source  approach  not  only 

) 

because  it  is  a  more  economical  approach  but  also  by  allowing  the  use  of 
a  larger  number  of  quadrilaterals,  the  results  are  more  accurate,  showing  a 
decisive  improvement  in  phases  and  moderate  improvement  in  magnitude  of  the 
pressure  points  compared  to  the  corresponding  experimental  values. 

As  for  the  "reflection  coefficient"  which  presents  the  effect  of  the 
image  due  to  the  presence  of  a  rigid  wall,  it  is  difficult  to  discern  any 
definite  trend  or  to  draw  any  conclusion.  The  commonly  used  approximation 
of  doubling  the  free  space  pressure  does  not  seem  to  be  valid  in  this  set  of 
calculations.  Conclusions  cannot  be  drawn  at  this  stage  until  additional 
calculations  are  performed  and  compared  with  experiments. 

It  would  seem  unlikely  that  a  reflection  factor  greatly  in  excess  of 
2.0  should  be  obtained  except  possibly  near  the  aft  extremity  of  the  skeg. 

In  any  event,  numerical  procedures  such  as  advanced  here  should  really  be 
compared  for  cases  in  which  an  independent  method  can  be  employed.  As  the 
Green's  function  for  a  spheroid  i  s  known  it  is  recommended  that  the  program 
be  employed  for  that  case  and  comparisons  made  with  the  outcome  of  the 
theory  for  a  spheroid  in  the  presence  of  a  single  rotating  source.  Only  in 
this  way  can  definitive  information  be  secured  regarding  the  maximum  panel 
size  permitted  to  achieve  a  requisite  accuracy. 
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APPENDIX 


EVALUATION  OF  A^z)  AND  A|n^(z)  FUNCTIONS 


n  -IzcosS 

A'n'  (z)  =  J  ©(n)  e  sin0  d© 

_ o _ 

a)  Birnbaixn  distribution 

A^^(z)  *  “j  J  cot  |  e  lzcos6sine  de  *  JQ(z)  -  iJj(z) 


A (n>l )  (z)  ■  ~  J  sin(n-l)0  sin0  e 


izcos© 


d0 


•  >  n  ”2 


[j-  <2> * j-(z)] 


n-2 


b)  "Roof-top"  distribution  (a  mean  lines) 


,(1).  .  pCOS  -izcos0  . 

A  (z)  «=  J  e  sin0  d© 


f"  (^c°Sgle~izC05S  sinB  dB 

J  _  2(1 -a) 

cos  (1 -2a) 


**  e 


■izji  .  i  f*  i2az  .  i2z"|  j 

|  z  2(1-5)  z2  L  J) 


(For  5  -  1.  A(,)(z)  .) 


A(l;>l)(z)  -  0 


c)  Sine  series  distribution 
,  (n)  *  ” 


A^(z)  ■  ~  J  sin  n©  sin©  e  '2C0S®  d0 
o 

.tfl  [j_  (I)+J.  W] 

2  L  n-1  n+1  J 
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(2) 


A^(z)  =  J*  .  @(n)  sin0  cosB  e  'zcos®  d © 

_ _ _ o _ 

a)  Birnbaun  distribution 

a}«W  '  Vz>]  '  iJl(z) 

a.<;;>,)(z)  (z)-j.  w] 

*  4  L  n-3  n+1  J 

b)  "Roof-top  distribution  (a  mean  lines) 


*<'>w 


ri  +  Jr  +  _LT[. 

'  z  z8  z2  (1-a)  L 


1 _  f  ( 1  -2a  -  1)  >2az 


(For  a  *=  1  ,  A 


<'><z)  (cos 


s  in 


*>•> 


c)  Sine  series  distribution 
kn+2 


A<"> (z) 


iiiX 


[■ 


J  (z) 

n-2 


J.  (z)l 

n+2 


It  is 
A^(-z)  and 


to  be  noted  that  the  values  for  negative  argument,  i.e., 
Ajn^(-z),  are  the  conjugates  of  the  values  given  above. 
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FIG.  3  DEFINITIONS  OF  ANGULAR  MEASURES 
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FIG.  4  RESOLUTION  OF  FORCES  AND  MOMENTS 
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FIG.  7.  Transducer  positions  on  skeg  and  transom,  projected 
[Fig. 3. 6  in  vertical  and  horizontal  plane  respectively.  Numbers 
in  circles  are  position  identifications.  Distances 

Ref. 

are  in  millimeters  on  model. 
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PROFILE  AT  i. 


CORRELATION  OF  MEASURED  AND  CALCULATED  PRESSURES 
(BOTH  WITH  SOURCE  DISTRIBUTION  AND  DOUBLET 
DISTRIBUTION)  .  PROFILE  AT  CENTERLINE. 
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FIG.  10  CORRELATION  OF  MEASURED  AND  CALCULATED, PRESSURES 
(BOTH  WITH  SOURCE  DISTRIBUTION  AND  DOUBLET 
DISTRIBUTION)  AT  SECTIONS  B-B  AND  D-D. 
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FIG.  11  CORRELATION  OF  MEASURED  AND  CALCULATED  PRESSURES 
(BOTH  WITH  SOURCE  DISTRIBUTION  AND  DOUBLET 
DISTRIBUTION)  AT  SECTIONS  E-E  AND  F-F. 


